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Abstract 



We consider a very simple Navier-Stokes model for compressible fluids in one space 
dimension. Neglecting the temperature equation and assuming both the pressure and the 
viscosity to be proportional to the density, we get 

dtipv) + d^ipv"^ + Aep)) = ed^{pd^v), 

dtp + d^{pv) = 0, 

where A, e are positive constants. We assume both the density field pit,-) > and the 
velocity field v{t, x) to be 1-periodic in space, with unit mean for the density. We denote 
by X(t,a) = a + ^(t,a) the "material" positions of the fluid parcels, so that ^ is one- 
periodic in a G M, dtX(t, a) = v(t, X(t, a)) and daX(t, a)p(t, X{t, a)) = 1. Then, we show 
that X(t, a) = a + ^(t, a) can be obtained as the unique limit, as h > goes to zero and 
nh goes to t, of the time-discrete approximation Xh,n{ci) = a + C,h,n{o,) obtained from the 
recursion formulae 

2) ^f^M = [a + (1 + hX)Ch,n-iia) + hZh,n-i{a)) + v^iV(a//i)]]» 

Zh,n{(^) = (1 - \h)Zh,n-iia') - A^^ft,„_i(a)), 

where Xh^Q, Z^^q are suitably initialized, is a 1-periodic function with average and 
variance 1, and ^ denotes the rearrangement operator in increasing order for functions 
depending on a G M. Since the term N{a/h) can be interpreted as a deterministic 
approximation of a random variable W with zero mean and unit variance, we see that, 
in some sense, the resolution of the Navier-Stokes system has been reduced, through the 
crucial intervention of the monotonic rearrangement operator, to the discretization of a 
trivial stochastic differential equation, namely 

= V2^W + {\^ + Z)dt, dZ = + Z)\dt. 



In addition, our result can be easily extended to a related Navier-Stokes-Poisson system. 
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1. Introduction 



Optimal transport theory has been succesfully used to treat many different types of 
parabohc equations as gradient flows of various "entropy functionals" for various "trans- 
portation metrics" , the canonical example being the regular scalar heat equation viewed 
by Jordan, Kinderlehrer and Otto [21j as the gradient flow of the Boltzmann entropy for 
the quadratic Monge-Kantorovich MK2 (frequently nammed Wasserstein) metric. There 
has been a huge literature on this subject in the last 20 years (study of the heat equation 
in a very general framework, porous-medium equations, thin-film flow equations, chemo- 
taxis models, etc., sample, [H |23l [2l [H]). Therefore, it is natural to try to 

treat the Navier-Stokes or the Euler equations of fluid mechanics in a similar way. Let 
us mention some attempts such as the introduction of a JKO scheme for the Euler equa- 
tion by Gangbo and Westdickenberg or the related treatment of the Navier-Stokes 
equation by Gigli and Mosconi [20] • However, we are probably still very far to be able 
to solve the initial value problem for the Euler or the Navier-Stokes equations by optimal 
transport tools. Nevertheless, there is more hope for fluid models in one space dimension. 
Indeed, in one dimension, optimal transport according to the MK2 metric is rather simple, 
as well known, since it can be entirely rephrased in terms of monotonic rearrangement 
and pseudo-inverse functions. Very simple models, such as the inviscid Burgers equa- 
tion (and more generally multidimensional scalar conservation laws) or one-dimensional 
pressure-less Euler equations, have been successfully written, using pseudo-inverse func- 
tions and monotonic rearrangements, directly as sub differential equations with nice con- 
traction properties in all spaces [21[ini[I2]- In the present paper, we show that similar 
ideas can be extended to some Navier-Stokes models for compressible fluids in one space 
dimension, leaving widely open the the much more challenging case of 2D and 3D Navier- 
Stokes equations. More precisely, let us consider the (over-simplified) fluid mechanics 
model of Navier-Stokes type, without energy equation, in one space variable: 



where p = p(t, x) > 0, v = v(t, x), respectively denote the density and the velocity of the 
fluid at time t > and point x, on the real line, and p, denote the pressure and the 
viscosity of the fluid that are supposed to be just given functions of the density. (The 
genuine Navier-Stokes equations would include an additional equation for the temperature 
field, on which both p and p would depend.) Our approach is restricted to the (rather 
unphysical) case: 



[As a matter of fact, if we neglect temperature effects, as we have done, it would be 
physically more consistent to assume the viscosity to be a constant and the pressure to 
be proportional to the density, but our method would not apply then.] The e parameter 
will be interpreted subsequently as a level of noise in the formulation we will propose for 
the Navier-Stokes system. In particular, the "pressure-less viscosity- less" case p = p = 
will correspond to a zero level of noise. The A parameter just scales the pressure with 
respect to the viscosity. 



(1.3) 



dtipv) + d^{pv'^ + p{p)) = d^{p{p)d^v), 

dtp + d^{pv) = 0, 



(1.4) 



M = ep, P{p) = p. 



Let us now write down the NS model (11.31) . under assumption (11. 4p . in "material" (or 
Lagrangian) coordinates. For that purpose, we define the material positions X{t,a) of 
the fiuid parcels which satisfy: 

(1.5) dtX{t,a) = v{t,X{t,a)), p{t, X{t, a))daX{t, a) = 1, X{t, a + 1) = X{t, a) + 1. 
Then, we find (after standard calculations with a multiple use of the chain rule) 

dtXit,a) = Vit,a)-dai ^^ ' J 

OaX(t, a) 

(1-6) X 

dtV{t,a) = -dai ' ' 



'daX{t,a) 
or, equivalently, 

X{t, a) = a + C,{t,a), 

(1.7) a) + daig-]^^) = Zit, a) + X^it, a) 

dtZ{t, a) = -\Z{t, a) - \^^{t, a) 

by substituting Z ioi V — X^. Our key observation is that this system looks like a very 
mild modification of 

X(t, a) = a + C,(t,a), 

^^■^^ mt,a) + d^{ ' )=0, 

daX{t,a) 

which is nothing but than the "material" version of the linear heat equation (written in 
"JKO style" pi]) 

(1.9) dtp + d,{pv)=0, v = -^, 

P 

through the same change of variables (11.51) . This suggests the following discrete scheme, 
based on the simulation of the heat equation by random walk, to approximate the Navier- 
Sokes equation written in material coordinates (II. 7p . 

Given a uniform time step h > and a uniform grid on M with mesh size 1/M. we 
approximate X{t,a), at each discrete time t = nh, n = 0, 1, 2, ■ ■ -, by a non- decreasing 
sequence k — )■ Xn± corresponding to each subinterval {k — 1)/M < a < k/M, for G Z. 
Because of spatial periodicity, it is convenient to introduce 

(1.10) ^{t,a) = X{t,a) - a, 

which is 1-periodic in a. Accordingly, we require Xn,k+M = 1 + Xn,k (which makes the 
scheme effective for computation) and introduce ^n,k = ^n,k ~ k/M. We first perform a 
predictor step 

^n+l,k = (,n+l,k + k/M 
(1-11) L+l,k = (1 + hX)^n,k + hZn,k + V2dlNn,k 



Zn+l,k — (1 — Xh)Zn,k — hX ^, 



n,ki 



where Nn^k is a sequence, M-periodic in /c, of M independent random numbers, with 
expectation and unit variance. So far, this can be seen as a discretization of the rather 
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trivial stochastic differential equation 



di={Z + \Cjdt + W 
(1.12) ^7 

where is a normalized white noise. However, we observe that there is no reason for 
k — )■ Xn+i^k to be increasing. So, we perform a corrector step just by sorting this sequence 
in increasing order, which is enough to define the next sequence k — )■ X„+i_fc. As discussed 
below, we even simplify the scheme by using a deterministic simulation of the white noise, 
for instance by substituting for the random number N^^k the value +1 if is even, and 
— 1 if k is odd, independently of n. This simplification is possible, because the corrector 
step generates a large amount of mixing. Let us point out that the sorting step is crucial 
to lift the rather uninteresting stochastic ordinary differential system fll.l2p to the level 
of the more intriguing Navier-Stokes equations! 

It is now very easy to define the natural continuous version of this scheme with respect 
to the space variable a, while keeping t discrete. To do that we use the concept of 
rearrangement of functions in increasing order, which is just the continuous version of the 
sorting operator for sequences of real numbers. However, we have to be a little bit careful 
to handle the required space periodicity. Let us consider a given function F G / + L^(T), 
where / denotes the identity map on M and L^(T) is the space of all locally square 
integrable 1-periodic functions on M. Then, its rearrangement in increasing order is the 
unique function E I + L^(T) such that 



:i-13) a„F«>0, fG{Y\a))da= f 

Jo Jo 



G{Y{a)da, 



for all 1-periodic continuous function G, G E G^{T). Then, the semi-discrete version of 
the scheme just reads (with obvious notation): 

Xn+i{a) = in+i{a) + a, 

4+1 («) = (1 + h\)Caia) + hZn{a) + V2dlW 

Zn+i{a) = (1 - \h)Zn{a) - /iA^^„(a), 

where is a normalized white noise. As a matter of fact, we rather use a deterministic 
approximation of the white noise term, namely 

(1.15) W = N{a/r), 

where iV is a fixed 1-periodic bounded function with zero mean and unit variance and 1/r 
is a sufficiently large integer. This prevents us to rely on any stochastic analysis. 
Let us point out, before stating our convergence result, that, due to its very definition, the 
semi-discrete scheme immediately becomes a fully discrete scheme, under the following 
conditions. We first approximate the initial data (Xo(a), Zo{a)) and the identity map a by 
piecewise constant functions in a, on a uniform grid, with constant values in each interval 
]{k — l)/M, k/M[, for k Eli. Observe that the uniformity of the grid is crucial, because the 
rearrangement operator in increasing order preserves piecewise constant functions only on 
a uniform grid. We also crucially need the "noise function" N{a/r) to be itself piecewise 
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constant on the same grid. Let us consider for instance the simplest noise function 
N{a) = -1 for < a < 1/2 and N{a) = 1 for 1/2 < a < 1. Since L = 1/r is supposed 
to be an integer, we see that N{a/r) will be piecewise constant on the grid of mesh 
width 1/M as soon as we choose M to be a multiple of 2L. If we choose, with economy, 
M = 2L, we see that N{a/r) will take value -1 on each interval ] {2k - 2) /M, {2k - 1) /M[ 
and +1 on each interval ]{2k — 1)/M,2k/M[, for G Z. Then the time-discrete scheme, 
without modification, coincides with the fully-discrete scheme discrete fll.lip where we 
substitute for the random number Nn^k the value -|-1 if /c is even, and —1 if /c is odd. Let us 
also observe that, by construction, the fully discrete scheme can be trivially coded (with 
the help of a fast sorting algorithm, taking into account periodicity), and the resulting 
complexity, to reach the solution at a finite fixed time t, is of order 0(/i~^r~^ log(l/r)) 
on a sequential machine. 
It is now time to state our main result. 



Theorem 1.1. Let us consider the time- discrete scheme {l.l^M-l^l . where N is the 1 



periodic binary function with value —1 for < a < 1/2 and +1 for 1/2 < a < 1, and 1/r 
an integer large enough so that r = 0{h). Let us consider initial conditions Xq = I + ^q, 
Zq, where I is the identity map and C,o and Zq are smooth 1 -periodic functions, with 
inf Xg(a) > 0. Let us linearly interpolate in time the {Xh.n, Zh,n) o,nd denote the result by 
{Xh, Zh). Then, as h ^ 0, the entire family converge to a solution {X, Z) of the material 
formulation of the Navier- Stokes system ( [i. Tp , under assumption 

Symbolically, to describe this limit, we suggest the following continuous formulation: 

X{t + dt, a) = [X{t, a) + {Z{t, a) + X{X{t, a) - a))dt + Wf 

(1-16) dZ 

^ + XZ = -X'{X-a), 

where W stands for a normalized white noise and " denotes the rearrangement operator 
in increasing order. 

Our method of proof will involve three main steps: 

i) The Navier-Stokes system, written in material coordinates, is first interpreted as a 
modified gradient flow with delay terms for a suitable convex functional, in a slightly 
more general situation. As a byproduct, we establish the global well-posedness of the 
system, in a general setting of non smooth initial conditions (including vacuum, i.e. zones 
where the density field vanishes). 

ii) In the special case (11.41) . we show that the previous system can be further reduced to 
a non-autonomous viscous scalar conservation laws, for which existence of global smooth 
solutions can be routinely established. 

iii) Finally, we observe that our time-discrete scheme is an adaptation of the "transport- 
collapse-kinetic" method [U [5l [191 El Ej to the discretization of a suitable coupled scalar 
parabolic equation (I3.37P and prove the desired convergence theorem. 

Let us conclude this introduction by mentioning that our approach can be adapted to 
some interesting cases when a self-consistent forcing term is added to the NS equations. 
This is in particular true for the Navier-Stokes-Poisson (NSP) system 
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dtp + d^{pv) = 0, f3dl^^ = p-l, 

where /3 is a constant, (j) = (j){t, x) can be interpreted ad a gravitational potential if 
/3 > 0, or as an electrostatic potential if /3 < 0, with a "neutralizing" background of unit 
mean. In material coordinates, this system reads 

(1.18) dlX{t, a) = da[-p{l/daX) + u{daX)dl,X{t, a)] + F{X{t, a) -a), 
where 

(1-19) F{y) = |. 

For this Navier-Stokes-Poisson system, under assumption ( II. 4p . we just have to modify 
the noisy differential system in a straightforward way 

d^={Z + X^)dt + V2^t W 
(1-20) dZ 

where F is given by fll.lQp . So we may immediately include the Navier-Stokes-Poisson 
case, just by adding a smooth forcing term, F{y) with bounded derivative in y and 
F(0) = 0, to the Navier-Stokes equations, as in fll.181) . without necessarily assuming that 
F has form fll.lQp . This is what we will do through the paper. 

2. Reduction to a gradient flow with delay 

In this section, we consider the Navier-Stokes model (written in material coordinates) 
fll.lSp . in the case when the viscosity coefficient /x and the pressure p are linked to each 
other through a given smooth strictly convex function r > — )■ 0{t), in the following 
way: 

(2.21) /i(l/r) = er^"(r) p(l/r) = -Ae ^'(r), Vr > 0. 



for some constants e, A > 0. We recover the case (11. 4p discussed in the introduction for 
the special choice: 

(2.22) e(r) = -logr, ^'(r) = -1/r, ^(r) = l/r^. 

Another interesting example is 9{t) = r(logr — 1), for which p = — Ae logr, and p = e 
(i.e. constant viscosity). Under assumption f l2.2ip . the NS equation written in material 
coordinates, with forcing F, (11.180 becomes 

(2.23) dlX{t, a) = Ae 4(^'(4X(t, a))) + e 9,(0" {daX{t, a))dlX{t, a)) + F{X{t, a) - a) 

which can be immediately blown up as a first order system: 

, , dtX{t, a) = V{t, a) + e da{e'{daX{t, a))) 

(2.24) 

^ ^ dtV{t,a) = Xeda{6'{daX{t,a))) + F{X{t,a)-a) 



or, equivalently 

X{t, a) = a + C,(t,a), 

(2.25) dt^{t, a) = Z{t, a) + X^{t, a) + e da{e'{daX{t, a))) 

dtZ{t, a) = -\Z{t, a) - A^^(t, a) + F(^(t, a)) 

just by introducing Z = V — X{X — a), i.e. V = Z + \{X — a). The resulting system 
is a mild modification of a gradient flow for a convex functional. More precisely, let 
us introduced the Hilbert space H = L^(T) and the closed convex cone K of all non- 
decreasing functions in I + H. 

Now let us introduce the Isc convex functional ^ E H ^ Q[^] e] — oo, +00], with value 
+00 if X = J + ,^ is not in K and 

(2.26) = / 0(1 + ^'{a))daE]- 00, +00], 

Jo 

otherwise. [Notice that whenever X = / + ^ G -ft', its derivative X'{a) = 1 + ^'{a) can 
be seen as a locally bounded nonnegative Borel measure on T. Then 0[^] can be more 
precisely defined as "a convex function of measures" by Legendre duality, as done, for 
example, in 



(2.27) e[e] = sup / {fia)d{a + ^a)) - 9*{f{a))da}, 

f Jt 

where / describes the set of continuous function on T and 9* is the Legendre-Fenchel 
transform of 9: 

9*{s) = snp{sT -9{t)). 

For example, in the case (12.221) we are mostly interested in, 9{t) = — logr (extended by 
+00 for r < 0), we find 9*{s) = — 1 — 9{—s) = — 1 — log(— s) and we can rewrite (12.271) as 

6[X] = sup / {— exp{g{a))d{a + ^{a)) + {g{a) + l)da}, 
g Jj 

for all continuous functions g on T, by setting /(a) = — exp{g{a)).] 

So, the first order system (12.251) can be viewed, in the Hilbert space H, as a sub-differential 
inclusion coupled to a linear ODE in Z, namely: 

-f + Z + A^G ede[^], 

(2.28) 

^ + XZ = -X'^ + F{0. 

By integrating out Z from the second equation, we may view the system as a gradient 
flow with delay terms (which are harmless since we have assumed F{y) to have bounded 
derivative in y and F{0) = 0, as in the Navier- Stokes- Poisson case (I1.19P ): 

(2.29) - ^ + Z{0)e-'' + l\-'^'-''\Fim,a) - \\{t'))dt' + G tdQii]- 

Such formulations are very useful to define global "generalized" solutions for a large 
set of initial conditions, with well-posedness in L^, just by using, in a routine way, the 
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general theory of maximal monotone operators |T3]. In addition, in this framework, the 
limit e — >^ is trouble-free. Indeed, in the sense of maximal monotone operator theory, the 
limit of eO[^] is just the indicator function of the cone with value whenever 

X = / + ^ belongs to K and +oo otherwise. Thus, we may conclude this section by 
asserting: 

Theorem 2.1. Let {^q^Zq) given in the Hilbert space H = L^(T), with Xq = / + 
in the closed convex cone K of all non- decreasing functions in I + H. Then, there is a 
unique global solution t — )■ (^(t), Z{t)) E H x H to l{2.28\) . depending continuously in H on 
both t and the initial condition (^(0) = C,o,Z{0) = Zq). This solution can be viewed as a 
"generalized" solution to the Navier-Stokes system U.18\) . written in material coordinates, 
with pressure and viscosity laws given by \2.21\) . and, in particular for OT^. In addition, 
as e — )■ 0, these solutions converge to those of: 

-^ + z + \ie diK[l + il 

(2.30) ^* 

^ + \Z = ~\'X + F{i). 
dt 

where is the indicator function of the cone K , with value whenever X = I + ^ 

belongs to K and +oo otherwise. 

Remark 1. Formulations of type ( 12.30p have already been introduced in the absence 
of viscosity for the pressure-less Euler equations and some of their variants. Let us 
quote the model of order-preserving vibrating strings, the Chaplygin gas, the pressure- 
less Euler- Poisson system and also multi-dimensional scalar conservation laws etc ... 
( [9] , [ID] , [22] , [H] , for example). In all of these earlier works, A has never been taken 
into account. It is unclear to us, whether or not A plays (or should play) a role in these 
zero viscosity limit equations. We leave that for further investigations. 

Remark 2. Subdifferential formulation (I2.28p . or (12.291) . leads to stability of solutions 
with respect to their initial condition not only in but also in all spaces. Indeed, let 
us consider, on a fixed finite time interval [0,T], two solutions {^,Z) and (^,2'). Using 
formulation (I2.29p . we get 

e\t) < c{||Z(0) - m\\l. + e{t) + fe{t')dt'], e{t) = \\m - m\h, 

Jo 

where c depends only on q, T, A and the Lipschitz constant of F, but not on e. This easily 
leads to the L'^ stability estimate, uniformly on e, 

(2.31) sup m,z)it) - iizmwL^ < cm,zm - iizmwL., 

te[o,T] 

with (another) c depending only on q, T, A, F, but not on e. 

3. Reduction to a scalar parabolic equation with mild coupling 

The sub-differential formulation fl2.28p of the Navier-Stokes equations in material co- 
ordinates (11.180 . that we have introduced in the previous section, is very powerful since 
it provides (and, more generally L*^) well-posedness for very general data (including 



vacuum). However, it is not clear that the corresponding solutions are smooth, for ap- 
propriate initial conditions. At least in the special case f l2.22p where 6{t) = — logr, this 
can be done by reducing fll.lSp to the scalar parabolic equation 

(3.32) dtu + {Z{t, u) + {x — u)X)dxU = ed^^u, 

where Z is coupled to m in a suitable way. For this purpose, let us introduce the inverse 
function u of X = I + 

(3.33) u{t,X{t,a)) = a. 

At the moment, let us assume X{t,a) = a + ^{t,a) to be smooth with daX > (which 
will be justified a posteriori in the case (11.41) ). We have: 

u{t,X{t,a)) = a, {d^u){t,X{t,a))daX{t,a) = 1, 
^' ' {dtu){t,X{t,a)) = -dMt,X{t,a))dtX{t,a). 

Let us now point-wise multiply the first equation of (I2.25P by 

{d.,u){t,X{t,a)) = l/daXit,a). 

Using (13.341) . we get 

~{dtu){t,X{t,a)) = [Z{t,u{t,X{t,a)) + X{X{t,a) - a)]{d^u){t,X{t,a)) 

+eda[e\l/id,u){t, X{t, a))]/daX{t, a), 
and therefore, by substituting x for X{t,a): 

—dtu{t, x) = [Z{t, u{t, x)) + A(x — u{t, x)]dxu{t, x) + edx[0'{l/dxu{t, x))]. 

Thus, we get for {u, Z) the following system: 

dtu{t, x) + [Z{t, u{t, x)) + A(x — u{t, x))]dxu{t, x) = —edx[9'{l/dxu{t, x))]. 
(3.35) „ 

dtZ{t,a) + XZ{t,a) = F{X~a)-X^{X{t,a)-a), u{t,X{t,a)) = a.. 

By integrating out Z, we get, more explicitly. 



(3.36) Z{t, a) = Z{0, a)e-^* + [ e-^(*-*')(F(X(t', a) - a) - X'^{X{t', a) - a))dt' . 

Jo 



Next, under assumption (I2.21|2.22l) . i.e. 9{t) = — logr, we find 

dtu(t, x) + [Z(t, u{t, x)) + X{x — u(t, x))]dxu(t, x) = ed^^u, 
^^'^'^'^ dtZ{t,a) + XZ{t,a) = F{X - a) - X'^{X{t,a) - a), u{t, X{t, a)) = a. 



Observe that the coupling of Z with u, through X is very mild, so that the theory for this 
system differs very little from the standard theory of viscous conservation laws, such as 

(3.38) dtu + dxifiu)) = edl,u. 

So, without entering in details, we claim that for uq given in the class C of all orientation- 
preserving diffeomorphisms of M/Z and Zq smooth 1-peridoic function, then there is a 
unique smooth solution {u{t,-) G C,,Z(t,-)) to the coupled scalar parabolic equation 
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stability with respect to initial conditions. In addition, thanks to the the subd- 
ifferential formulation f l2.28p . we get an stability result for every pair {u, Z) and {u, Z) 
of solutions to the coupled scalar parabolic equation f l3.37p . on every fixed finite time 
interval [0, T]: 

(3.39) sup (|Kt)-^(t)|U, + ||Z(t)-Z(t)|UJ <c(|K0)-i2(0)|U, + p(0)-Z(0)|UJ 

te[o,T] 

with c depending only on T, A, F, but not on e. This immediately follows from the L'^ 
stability estimate f l2.3ip for the subdifferential formulation fl2.28p . in the special case 
q = l. 

4. Approximation by a "transport-collapse-kinetic" method with noise 

The transport-collapse-kinetic method is a time-discrete scheme for hyperbolic and 
viscous conservation laws introduced and studied in [U [191 El [8] . It was suggested to the 
author by an earlier algorithm proposed by A. Chorin [16] for reaction-diffusion equations. 
This scheme turns out to be very well suited for the formulations f l3.37p and f l2.28p we 
have obtained in the previous sections for the Navier-Stokes model fll.l8p . 
Let us now describe the time-discrete scheme: 

Given a uniform time step h > 0, we denote by (X/j „(a) = a + ^h,n{(^), Zh,n{ci)) the 
approximate solution at time t = nh, for = 0, 1, 2, ■ ■ defined in two steps as follows. 
We first define Zh^^+i and Xh,n+i = a + 4,n+i(a) by 



4n+i(«) = (1 + hX)^hAa) + hZh,n{a)) + V2dlN{a/r) 

Zh,n+l{o) = (1 - Xh)Zh,n{a) + h{F{S,h,n{o)) - X^ih,n{o)), 



where is a 1-periodic function with average and variance 1, and 1/r is a large integer 
that may depend on > 0. [Let us recall that this first step has been designed to 
approximate the noisy ordinary differential system (11.120 . namely: 



di = {Z + XC)dt + W 
^'■''^ f + AZ = F(0 - \% 

where is a normalized white noise, with a deterministic approximation to W ] 
Next, we crucially introduce a rearrangement step. Namely, we rearrange Xh^n+i{o) = 
a + ^/i,n+i(a) in increasing order with respect to a, and obtain Xh,n+i{,o) = a + S,h,n+i{,o). 
In other words, we set 

(4-42) X,,„+i = 

where we denote by y — the rearrangement in increasing order, defined by (11.130 . for 
functions differing from the identity map / by a L'' 1-periodic function (i.e. belonging to 
/ + L'^iT)). Of course, the resulting scheme contains the scheme (I1.14p discussed in the 
introduction, up to the addition of the (harmless) term F, that we suppose, as previously 
mentioned, smooth with bounded derivative. Notice that the scheme is almost translation 
invariant in the variable a. Indeed, the rearrangement step is clearly translation invariant, 
and, in the predictor step, the only non-autonomous term is the "noise" term which is 
proportional to N{a/r). Thus the scheme is invariant by any translation a — )■ a + k/r, 
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with /c G Z. This observation will be very important for our analysis. We also observe 
that Zft, is entirely "slaved" by and depends linearly on it through (14 .40 p . so we are only 
concerned about monitoring C,h- (Of course we use that F is smooth and Lipschitz with 
-F(O) = 0.) For the subsequent analysis, it is also simpler to restrict ourself to the case 
when is the binary function with value —1 for < a < 1/2 and +1 for 1/2 < a < 1. 
Indeed, in that case, remarkably enough, {^q, Zq) = (0, 0), is a fixed point of the scheme, 
as F = 0. [This can be easily checked graphically: we see the exact compensation between 
the addition of the "noise" term \/2ehN{a/r) and the rearrangement step.] 

Let us now prove the following convergence result (from which our main result Theorem 
II. II immediately follows as the particular case when the force term F is absent): 



Theorem 4.1. Let us consider the time- discrete ^^J^JUj^J^J^ , where N is the 1-periodic 



binary function with value —1 for < a < 1/2 and +1 for 1/2 < a < 1 and 1/r an 
integer large enough so that r = 0{h). Let us consider initial conditions Xq = I + C,o, 
Zq, where I is the identity map and o;^^^ Zq are smooth 1-periodic functions, with 
infXQ(a) > 0. Let us linearly interpolate in time the output of the semi-discrete scheme 
{Xh,n, Zh,n) CLnd denote the result by {Xh, Zh)- Then, as /i — )■ 0, the entire family converge 
to (X, Z) solution of the material formulation of the Navier-Stokes system U.18\) . under 
assumption (T^. 

Proof. 

Sup-norm and Lipschitz estimates. The rearrangement operator is a non-expansive map 
on all spaces / + L'^{T), in particular for q = oo. So we can compare in sup-norm any 
fixed solution {^h,n, Zfi^n){(^) of the time-discrete scheme to: 

i) (^h.n, ^/i,n)(« + kr), its space translation by kr, for any integer k E Z, which is also 
solution of the scheme (due to translation invariance, as already mentioned); 

ii) the fixed point solution already {C,,Z) = (0,0), obtained in the case F = 0, Let us 
perform the first comparison. We easily get: 

(4 43) ~ (eM,^M)(- + MI|oo < (l + M"ll(eo,^o) - (eo,^o)(- + fcr)|U 

< \k\rLip{^Q, Zq) exp{nhc), 

where Lip denotes Lipschitz constants and c is a constant depending only on A and the 
Lipschitz constant of F. [Indeed, we first get, for the predictor step, 

||(^/i,n+l, Zh,n+l) — {C,h,n+1, Zh,n+l){- + ^'")||oo 

<{l + ch)m 

,n+l) Zfi^n+l) {^h,n+ly ■2'/i,n+l) (' + kr) \\ooy 

where e and the "noise" term do not play any role, since the translation is an integer 
multiple of r. Next, using the contraction property of the rearrangement operator, we get 

which is enough to get (14.431 ] As a matter of fact, we can improve estimate (14.431) using 
the fact that Xh,n{cL) = a-\-^h,n{ci) is increasing in a. Indeed, we get, for all u G [kr — r, kr], 

ih,n{a + uj)- ih,n{a) = Xh,n{a + u) - Xh,n{a) < Xh,n{a + kr) - Xh,n{a). 

Thus, 

\\^h,n- (^/i,n(" +w)l|oo < \k\rLip{^o, Zo) exp{nhc) < {\u\ + r)Lip{^o, Zq) exp{nhc). 
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Since Z^^n is slaved by we have obtained 

(4.44) \\{ih,n,Zh,n) - {ih,n,Zh,n){- + (^)\\oo < (|w| + r)Lip(^o, ^"0) exp(n/ic) , 

for all a; G R, where c is a constant depending only on A and the Lipschitz constant of 
F and not on e. This is our first key estimate for the scheme. Next, let us compare the 
solution to the fixed point of the scheme (^ = 0, Z = 0), already mentioned and obtained 
for F = 0. For this fixed point, the predictor step provides the values {\/2ehN{a/r), 0) 
Thus, we get 



,n+li Zfi^n-\-i 

)-(o,o)iioo< ii(a oo 

(because the rearrangement operator is non-expansive) 

<(i + c/i)ii(a ,ni ^h,n) I |oo) 

(by definition of the predictor step, using that F is Lipschitz and F{0) = 0), where c is a 
constant depending only on A and the Lipschitz constant of F. So, we have obtained our 
second important bound 

(4.45) ||(a,n,^M)lloo < 11(^0, ^o)||ocexp(n/ic), 

where c depends only on A and the Lipschitz constant of F, and not on e. 

Consistency estimates. From now on, we fix a time interval [0,T] and denote by c any 
constant depending only on the data T, Xq = / + ^o, Zo, N, F,e, X. Let G be a smooth 
1-periodic test function. Let us fix a time step n. Since G is 1-periodic, because of the 
rearrangement step, we have 

(4.46) / GiXh,n+iia))da= [ G{Xh,n+i{a))da. 



Thus, using the definition of scheme (l4.4Up . we get 



G{Xh,n+i{a))da 

(4.47) 

G{Xh,n{a) + /i[Aa,n(a) + ^h,„(a)] + V2dlN{a/r))da. 



By Taylor expansion in h, we get for the right-hand side of f l4.47p . 

{GiXh,ni(^)) + G"(X^,„(a))/i(Aa,n(a) + Z^,„(a)) + G'iXh,nia))V2dlN{a/r) 



+^G"'(X,,„(a))[/i(Aa,n(a) + Z,,„(a)) + V2dlN{a/r)Y}da + ©(/i^/^), 

using the sup-norm estimate fl4.45p . We want to simplify this expression, in order to 
obtain 

»i 

{G'(X^,„(a)) + G"(X;,,„(a))/i(Aa,n(a) + ^/^,n(a)) 



+G" {Xh,n{a))eh}da + Oi^/^), 
There are two terms to deal with: 

Ii= [ G'{Xh^n{a))V2dlN{a/r)da 
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and 
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h = - G"(X,,„(a))[/i(Aa,n(a) + Z,,„(a)) + V2ehN{a/r)Yda, 
^ Jo 

which both involve v^- Let us write N(a) = M'(a) where M is a Lipschitz 1-periodic 
function Q, so that N{a/r) = rda[M{a/r)]. Thus, integrating Ji by part in a, we get 
Ji = -J3 + h 



h = - G''iXh,n{a))daXh,n{(^)V2dirM{a/r)da 
Jo 

and the boundary term 

h = rv^[G"(X;,,„(l))M(l/r) - G'(X;,,„(0))M(0)]. 

Since 1/r is assumed to be an integer, by 1-periodicity of G, Xh,n — I and M, we first 
find Is = 0. Then, using that X/j^„(a) is non- decreasing in a, we get 

I/3I < cr^/h / daXh,n{a)da = cr\/^(l + / daih,n{a)da), 
Jo Jo 

where c denotes a generic constant depending on the data and the test function G. Thanks 

to the sup norm estimate (14.451) . we finally get |/i| < cry/h. Thus Ji can be absorbed in 

the error term 0(/i^/^), just by assuming r = 0{h). Let us now consider I2 = I5 + Iq + Ij, 

where ^ 

h = l [ G"(X,,„(a))[/i(Aa,n(a) + Z,,„(a))]2c?a 
^ Jo 

1 

G"'(X^,„(a))/i(Aa,n(a) + Zh,nia))V2dlN{a/r)da 

1 

thG'' {Xh,n{a))N^{a/r)da 

Since we know by ( I4.45P that (^/i, Z^j is bounded in we get + Iq = 0{h^). The last 
term, J7 is just 

-1 

e/iG" {Xh^ri{,o))da, 



since N^{a) = 1. Thus, assuming r = 0{h), we have obtained the key consistency 
property of the scheme 

]- [\G{Xj,^n+i{a)) - G{Xh,n{a))]da = 

(4.48) _ 

' {G"(X;,,„(a))(Aa,n(a) + Z^Aa)) + eG" {Xh,n{a))}da + 0{Vh), 



for all smooth 1-periodic function G. We can write this relation in terms of the probability 
measure ph,n{dx) and the real- valued measure Qh,n{dx) defined by 

/ 9{x)ph,n{dx) = I g{XhAa))da, 

(4.49) 

g{x)Qh,n{dx) = I 9{Xh,n{a)){Kh,n{o) + ZhAa))da, \fg e C°(T). 
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Notice that both are uniformly bounded in {h, n) as Borel measures, since X/j „ and X^^^ 
are uniformly bounded, as nh < T, for all finite T. We obtain, as a substitute for (14.481) . 

(4.50) P^,n+l - PH^n ^ g^Q^^^ ^ ^ ^^^^^ 

in the sense of distributions on T. 

Compactness. We first perform a linear interpolation in time of the measures {ph,ni Qh,n) 
which generates a continuous piecewise linear function of time t — {ph{t), Qh{t)) valued in 
the dual space C°(T)' ( i.e. the space of all bounded Borel measures on T, by Riesz iden- 
tification theorem). For each fixed T > 0, this family is bounded in Lip([0, T], C^(T)'), 
because of fl4.50p . and therefore compact in C°([0, T], C°(T)'). Let us consider a conver- 
gent subsequence, for a sequence of time steps hm going to zero, and denote the limit 

(P,Q). 

Since Zh is "slaved" by C,h and depends "hnearly" on it through (14.401) . we are just 
concerned by the compactness of C,h- Because of the sup norm and Lipschitz estimates 
f l4.45f4.44p . we already know that (^/i(t)) is valued in a fixed compact set of L^(T). Next, 
we use in a crucial way the following lemma which comes from rearrangement theory (or, 
alternately, from "optimal transport theory") [71 [23]: 

Proposition 4.2. Let {^h) a family of functions in H = L^(T) such that = I + C,h is 
a non- decreasing function. Then {C,h) converges to ^ in H if and only if 

(4.51) [ G{Xh{a))da^ [ G{X{a))da, 

Jo Jo 

where X = I + ^, for all smooth 1-periodic unction G. 

Using this proposition, we deduce from fl4.48p that {C,h) is uniformly equi-continuous 
from [0,T] to H, for all < T < +oo. We conclude, by the Arzela-Ascoli theorem that 
the family {C,h,Zh) is relatively compact in C°(]R_|_, if), remembering that Zh is linearly 
slaved by C^h- 

Consistency. We now have enough compactness to pass to the limit: up to the extraction 
of a subsequence of h 0, we have at least a limit {^,Z) in G^{M.,H) for the family 
{^h, Zfi). Because of the sup-norm and Lipschitz estimates fl4.45|4.44l) . we also know that 
{^{t, a), Z{t, a)) are bounded and Lipschitz continuous in a, uniformly in t for all bounded 
time interval t G [0,T]. Concerning Zh, passing to the limit in the second equation of 
(I4.40p is straightforward (using that F is smooth and Lipschitz). We get: 

(4.52) ^ + = - 

Next, from ( I4.48p . we deduce that for all smooth 1-periodic function G 

d 

— / G{a + ^{t,a))da = 

(4.53) 

{G'{a + ^(t, a))(A^(t, a) + Z{t, a)) + eG" (a + ^(t, a))}da. 
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Uniqueness of the limit. Let us consider any limit Z) of the scheme, on a finite interval 
of time [0,T]. We know that it belongs to L°°{[0,T], Lip{T)). Thus, for each t, X{t,a) = 
a + ^{t,a) has an inverse function u{t,x) which is increasing and satisfies u{t,x + 1) = 
1 + u{t,x). In addition, there is a constant a > such that dxu{t,x) > a. (However, at 
this stage, u{t, x) may have jumps in x). Its derivative p(t, x) = d^uit, x) > a is, for each 
t, a probability measure on T, bounded away from zero. We may also introduce, for each 
time t, the real valued bounded (uniformly in t) measure Q{t,x) on T, defined by 

(4.54) f g{x)Q{t,dx)= [ g{X{t,a)){X^{t,a) + Z{t,a))da, \/g e C%T), 
Jo Jo 

which is absolutely continuous with respect to p(t, x) with a bounded density (since X 
and Z are bounded). Thus, (p, Q) belongs to L°°{[0, T], C°(T)'), where C°(T)' is the space 
of all bounded (Borel) measures on T. The consistency relation fl4.53p . written in terms 
of p and Q, exactly means 

(4.55) dtp + d^Q = ed^^p, 

in the sense of distribution. Their Fourier coefficients Fp{t,k), FQ{t,k) are uniformly 
bounded in t G [0, T] and k E Z, and satisfy 

(4.56) {dt + 47Thk^)Fp{t, k) = 2t7TkFQ{t, k). 

This implies that p belongs to L°°([0, T], if'^(T)) (where if*(T) denotes the standard 
Sobolev space of 1-periodic function with derivatives in up to order s, for s G M) 
for all s < 1/2, and, therefore, by Sobolev embedding theorem to all L°°([0, T], L'^(T)) 
for 1 < g < +00. Since the density of Q with respect to p is bounded, we deduce 
that Q G L°°([0, T], L'^(T)) for all 1 < g < +oo. Since p = dxU, we also deduce u G 
L^{[0,T],H-'{r)) for all s < 3/2, and, thus u G L°°([0, T], ^"■'^(T)) for all a < 1. Thus, 
from the definition of Q fl4.54p . we can now write 

(4.57) p{t,x) = dxu{t,x), Q{t,x) = {{x - u{t,x))X + Z{t,u{t,x))dxu{t,x) 

and deduce that Q, just like p, belongs to L~([0, T], iJ"(T)) for all s < 1/2. [For that, 
instead of using Fourier coefficients, we estimate translations of Q in L^. Using that u is 
bounded and Z Lipschitz, we get 

»i 

\Q{t,x + u)~Q{x)\'^dx < {h + hy 

Jo 

where ^ 

-^1 = / \p{t,x + uo) — p{x)\^dx 
Jo 

\u{t,x + u) — u{x)Y'dx I p(t,xydx. 
'0 Jo 

Since p belongs to L°°([0, T], //"*(¥)) for all s < 1/2, we immediately get that 

h < Cs|wp", for all s < 1/2. Next, we use that u belongs to L°°([0, T], if"(T)) for all 

s < 3/2. This implies by Sobolev embedding theorem, that u and Z{u) certainly belong 

to L°°([0, T], W^'\T)), while p belongs to L°°([0, T], L\T)). This implies that h < c\uj\\ 

So we see that / < Cs|a;p'' for all s < 1/2, which, indeed, means that Q belongs to 
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L'^{[0,T], H^{T)) for all s < 1/2.] So we can now bootstrap the regularity of p and Q, 
using fli36|) . We first get that p, belongs to L^{[0,T], H'(T)) for all s < 3/2 and u to 
the this space for all s < 5/2. By Sobolev embedding theorem, this means that u also 
belong to L°°([0, T], ^^''^(T)) for all a < 1. We already know that d^u{t,x) > a > 0. 
Since X{t, a) is the inverse function of u, it follows that ^{t, a) = X{t, a) — a also belongs 
to L°°{[0,T],C^''^{T)) for all a < 1. This is also true for Z, which is slaved by ^ through 
fl4.52p . At this point, the boostrap argument becomes pure routine and we conclude that 
u,^,Z are smooth, just as the initial conditions. From (14. 55114. 57f4.52p . we see they are 
the unique solution in classical sense of the system 

u(t,X(t,a)) = a, X{t,a) = a + ^(t,a), 

(4.58) dtu{t, x) + [Z{t, u{t, x)) + A(x - u{t, x))]dxu{t, x) = edl^u, 

dtZ{t, a) + \Z{t, a) = F{^{t, a)) - \^^{t, a). 

This concludes the proof of our theorem since this system is just the formulation (13.371) 
of the Navier-Stokes system (11.181) . 
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